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Abstract 

We develop an exact generalized Bogoliubov transformation for the spin 3/2 Hubbard model 
with large anti-Hunds rule coupling near half filling. Since our transformation is unitary, we can 
thereafter employ standard approximate mean field theory methods in the full Hilbert space to 
analyze the doped Mott insulator, in contrast to a conventional approach based on truncated 
Hilbert spaces complemented with hard core constraints. The ground state at exactly half filling 
is an insulating (Mott) singlet, and according to our analysis a non-Fermi liquid order parameter 
A usually associated with extended s-wave superconductivity, will appear self-consistently as soon 
as a finite density n holes are introduced. The non-Fermi liquid behavior is a consequence of the 
nonlinear nature of the unitary transformation mapping the Mott singlet state to a Fock vacuum 
which introduces anomalous terms such as An in the effective Lagrangian. Our analysis uses an 
approach that generalizes readily to multi-band Hubbard models and could provide a mechanism 
whereby a non- Fermi liquid order parameter proportional to density is developed in Mott insulators 
with locally entangled ground states. For more complicated systems, such an order parameter could 
coexist naturally with a variety of other order parameters. 
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I. INTRODUCTION 

A Mott transition is expected to occur when the overlap between atomic orbitals in an 
insulator becomes large enough for the hopping energy to overcome the energy associated 
with charge fluctuations. To understand this transition has turned out to be a very difficult 
problem that has been attacked by a variety of means, with dynamic mean field theory 
being an important recent contribution |l|, |2j, |3|, y, |5|, la, [7|. In the present work we will 
develop an alternative approach, more in the spirit of BCS theory, which is appropriate for 
understanding doped Mott insulators where the parent state has an even number of electrons, 
and thus integer spin, at each state. Such models have been considered for instance in two- 
band Hubbard models in the context of ruthenate alloys 8| , and multi-band Hubbard models 
in the theory of doped Cqq J3 IM la l2l LLi| ■ We will present evidence for non- Fermi liquid 
behavior in these systems even for weak doping. We shall focus on the simplest system to 

namely the s pin 3/2 Hubbard model on a square lattice with 



ilalnilii, 



which our conclusions app 

anti-Hunds rule couplings [$, |£J 110", llll 112 Hi 

We define the local particle number by n r = J2 S c l r c s r an d the spin by S r = 
^2 s cl r S SjS ,c s/r , where — | < s < § and S s , s , are the generators of spin 3/2 rotations. Fur- 
thermore, we define the operator Pl m (r) which creates an / = 2, l z = m state with two 
fermions |l2j, P^ m { r ) = S^sd l^lf fi 2, rn )c [ a ^ r c f3r . We also define the SU(2) invariant 
P 2 = X] m ^2m( r )-P2,m( r )- The Hamiltonian containing the maximal number of onsite terms 
permitted by symmetry is 

H = -tJ2 4, r c s , r+d + X)(tMnr - 2) + JP?) , (1) 

r,S r 

where an arbitrary chemical potential is absorbed in U . Another term which can be consid- 
ered is of course S^, but for spin 3/2 this obeys P r 2 = ^S 2 — \n + |n 2 , so there is no loss 
of generality in the onsite term of our Hamiltonian aside from ignoring the possibility of an 
arbitrary term proportional to n 3 . We will consider the case J > 2U ^> t, thus the singlet 
state is heavily favored near n = 2. 

The single site spectrum consists of sixteen states: an empty site, four equivalent spin 
3/2 singly charged states, a singlet and five spin 2 doubly charged sites, four spin 3/2 charge 
three states and a charge four singlet. If E g (n) is the atomic ground state for n particles we 
find E g (0) = E g (2) = 0, E g (l) = -U, E g (3) = 5 J + 3U, E g (A) = 10 J + W and the quintet 



state has energy 2.1 . The lowest energy states obey E g (n+ 1) +E g (n— 1) — 2E g {n) > U 3> 
so there is no tendency for superconducting pair formation from any of the local interactions. 

With our choice of parameters, standard arguments imply that the ground state for n = 2 
and small t should be a spin singlet. For small doping near zero filling, the ground state 
will most likely be a normal Fermi liquid. This statement may be violated if U and J are 
sufficiently large, and a spin-symmetry breaking state may appear according to the Stoner 
criterion. Whether or not this happens depends on the density of states, which at least 
in the case of two dimensions remains finite even down to n = 0. Numerical simulations 
indicate the tendency toward spin ordering is grossly exaggerated in mean field theory [14]. 
In any case a deeper discussion of this point is not within the scope of the present article, 
which is to investigate the analog of a Fermi liquid near half filling. 

The Mott singlet | $ s (r) ) = A' r | ) at site r is created from the vacuum by the operator 
A\ given by 

At -^E e ^ ,c !A- (2) 

* s 

A natural first attempt to understand the system for small hole doping is to try the same 
method as for the nearly half filled spin 1/2 Hubbard model, i.e., to make a particle- hole 
transformation c\ r = c_ sr , where c\ r and c s r are new local fermionic creation and destruc- 
tion operators. In contrast to the case of the filled spin 1/2 Hubbard model however, this 
canonical transformation fails a number of criteria if the Mott singlet is to act as a vacuum 
for the new operators. In particular, we see that 

c., r | *.(0 ) * (3) 

( $ s (r) \c Str c{ r \ *.(r) ) + 1 • (4) 

The first inequality is perhaps not so severe; after all the operator c sr , being fundamen- 
tally a fermion creation operator creates a state with three fermions on site r, and we could 
argue that we could ignore this problem by suitably projecting the states on those with two 
or less particles. This is in fact the approach usually taken in attempts to perturbatively 
construct a new vacuum at n = 2. The second inequality is much worse; it is a consequence 
of the singlet being entangled, i. e. it cannot be written as a product state in any one-particle 
basis. As a result, the putative creation operator does not generate a normalized state from 
the vacuum. The entanglement property implies that the destruction operator on the ere- 



ated state does not recreate the ground state — not even in the two particle space since it 
projects into the 5 = 2 two-particle states. 

II. CANONICAL TRANSFORMATION TO THE MOTT SINGLET VACUUM 

We now show how to systematically construct creation and annihilation operators that 
have the correct local properties. We seek a canonical transformation that fulfills the fol- 
lowing criteria: (a) it maps | $ s (r) ) to | ), and (b) it maps the singly charged states to 
themselves. Due to our choice of interaction parameters, where J > 2U, we also expect that 
the state composed of two holes will play the same role in the hole doped Mott insulator 
near n = 2, as the doubly charged singlet does for small filling. The canonical transforma- 
tion that we desire therefore has the property that it interchanges the Mott singlet and the 
vacuum, leaving all other states invariant. The 5 = 2 doubly charged states will have the 
same charge as the Mott singlet. However, due to the constraints imposed by a canonical 
transformation utilizing only the spin 3/2 fermion operators, these doubly charged states 
must be obtained by two applications of the new creation operator, while the singlet cre- 
ated by another double combination of these creation operators forms a state with relative 
charge minus two. These considerations force the canonical transformation to be charge 
non-conserving. 



The canonical transformation that accomplishes t 



lis and similar mappings can be sys- 
However, we can obtain it without 



tematically obtained through the method in Ref. [15|. 
much formalism through the following argument. Our desired operator is "almost" Aj. 
The problem is that this operator generates unwanted side effects in the n = 1 and n = 2 
particle subspaces by mapping these to new states with n = 3 and n — 4. We can get 
rid of these unwanted overlaps by using a projection operator, and therefore define Q\ 
by Ql = A t r (l - n)(l - f). It is straightforward to check that Q\\ ) = | $ s (r) ), that 
Q r \ $ a (r) ) = | ), and that Q r and Ql annihilate all other states. A canonical transforma- 
tion that rotates the states | ) and | $ s (r) ) at each site r into each other without affecting 
the other states is provided by the unitary operator 

r 

U r (4> r ) = e iGr ^' r) = jM+faQr) ( 5 ) 



where <p r = (f) r e %Xr (<f) r and Xr are real). On | ) and | $ s (r) ) the transformation becomes 

U r {<j> r )\o,9 = cos <f> r + i sin (f) r {e iXr Ql + e- iXr Q r ) , (6) 

whereas it is unity on all other states. Choosing <p r = n/2, for all r, we obtain the canonical 
transformation that fulfills our criteria, i.e. it interchanges the empty state and the Mott 
state at each site without affecting the other states. 

Applying the unitary transformation Eq. El with r = tt/2, the true vacuum state, | ), 
is mapped onto the Mott insulator at half-filling | $ s ). The phase factors e %Xr enter this 
state only as an overall phase ^ r Xr, and can be neglected. In general, however, it is obvious 
from Eq. EJthat the unitary transformation gives a state where the phase factors enter in a 
non-trivial way. In particular, this is the case for the slightly doped Mott insulator, which 
we will be interested in below. This will be mapped onto a state near the true vacuum state, 
which can then be analyzed with standard methods. 

Note that the phase factors, e* Xr , are crucial to retain local gauge symmetry in the same 
way as the complex phases introduced into the Bogoliubov transformations are necessary to 
restore gauge invariance in BCS theory. 

III. A VARIATIONAL ANSATZ 

We now turn to a systematic variational analysis of the slightly doped Mott insulator using 
the canonical transformation in Eq. and use computer algebra to handle the complicated 
fermion polynomials that occur. In analogy with ordinary Fermi liquid theory, as well as the 
BCS theory of superconductivity, we will then search for a variational state with particle 
number given by n = 2 — 5 that is obtained from the vacuum by a canonical transformation 
e iQu depending on a set of parameters u. We define the functions E{u) and N(u) by 

E(u) = (0 \e- iGu He iG "\ 0) (7) 

N(u) = (0 \e- i5u Ne iGu \ 0) . (8) 

The values {u} which minimize E(u) define our variational ground state e lGu \ ) with 
particle number N(u). We have seen that for n = 2, the transformation Q u is simply 
G(|e iXr ) = G (xr) given by Eq. 0. We therefore expect that near the Mott insulator, the 
relevant transformation will be given by a further transformation close to the identity. We 
therefore make the ansatz e lGu = e 4G °( Xr V e " = U U'. 



Note that since we can continuously rotate the Mott state at n = 2 to the true vacuum 
by letting go from it/2 to zero, we can generate a Mott singlet on a site either by having 
U' r — 1 and <p r = 7r/2, or by having U' r = C/ r (7r/2) and r = 0. In general, we can make a 
coherent superposition of empty and doubly occupied singlet sites, both by letting <fi vary 
and by adding an onsite s-wave order parameter. As could be expected, this indeterminacy 
leads to a numerical instability in the variational equations which we resolve by simply 
taking <p r . = tt/2 for all r, and not further exploiting these variational parameters. 

In order to construct an ansatz for e l<3 ' u we first work out e~ tG °^ Xr ^ He tG °^ Xr \ This operator 
is obtained by replacing each occurence of the fermion operator cj s by e^ lGQ ^ Xr ' c\ s e %G °^ Xr ' 
and similarly for c sr . This expression is complicated, but it can nonetheless be worked out 
exactly in terms of polynomials of c sr and c\ r1 since the fermion algebra at a site is closed. 
The exact expression, written here for reference only, is given by 

4, r -► < r ((Ate- 2 ** (1 - n) - e 2ix A) + [^ + n + ^)) + 

(_l)( S +|) e 2,x 2 -| (-1 + e -2*x A t + (■£■ + n + ■£■)) c _ sr , ( 9 ) 

where the subscripts are dropped on the right hand side. The notation : O : is used to 
indicate a normal ordered operator, i.e. strings of fermion operators where all creation 
operators are anticommuted to the left and annihilation operators to the right taking only 
into account the sign of the permutation. In this case, : n 2 := n 2 — n and : S 2 := S 2 — 15n/4. 
The onsite interaction is zero in the vacuum and two particle singlet subspace. Since these 
are the only two states affected by the canonical transformation, this interaction remains 
invariant, while the chemical potential transforms according to 

rw 2-(n-*£i-^ + 3£) . (10) 

Anticipating a mean field calculation under the assumption of no spontaneously bro- 
ken global symmetries, we do a Wick decomposition of the onsite term, and calculate the 
expectation value according to 

(Un(n-2) + JP 2 ) = -hU + n 2 f^l + ^pj+2A 2 U , (11) 

where a hat indicates the expectation value of an operator composed of ordinary fermion 
operators evaluated in the state U'\ ) near the physical vacuum. Similarly, the expectation 



values for the density and s-wave order parameter AV become exactly 

(At) = ^£e--Ae^ (l + £!I - \ + £) - A' (|A|» -§ + £') (12) 

n \ I _ n \ n A T 



We now derive a similar expansion of the hopping operator. In this case the expressions 
become a terrible mess — far to complicated to write down in their entirety. It is however 
possible to construct a systematic expansion in the number of fermion operators, which is 
appropriate for small doping. To this respect, we take the entire expression, and rewrite it 
exactly as a sum of normal ordered terms. We then truncate this expression at fourth order 
in fermion operators and keep expectation values of on-site and nearest neighbor s-wave 
pairing amplitudes, ordinary hopping and density operators. Defining 

At - v t_ll_!_!_ c t c t ( U ) 

h r,r> = E s c{ r C s y , (15) 

we find that the hopping operator, truncated to fourth order in fermion operators, becomes 

( h r , r , ) = -f h (e 2 ^ A ry + e- 2l Xr'A* ry ) + \ e 2i ^-x r ,) (i _ fi) h r ,, r - 

t (e~ 2iXr A* r - e 2iXr A r ] e- 2ix r>A* rr , + 1 (V 2 ^' A;, - e 2ix -' A r >) e 2i ^k r y + ... (16) 

Decomposing the s-wave order parameters as, 

A r y = e~ iv ^'A r y (17) 

A r = e-^A r , 

where A r y is real, we obtain, 

( h ry ) = -te i{Xr ~ x ^ [\n(e i{Xr+x ^'^y)^ ry + h. c .) - 1(1 - n)e i( - Xr - x ^h r y (18) 

+ \(e 2i ^" r ^A r - h.c.)(e i( - Xr+x r'^r,r')A r y + h.c.)] + ... . 



Written in this way, the invariance under local gauge transformations is manifest (of course 
provided that the original hopping term is supplemented with the usual electromagnetic 
phase factor exp(i f dr" • A(f"))). 

Note that the combinations Xr + Xr' ~ Vr,r' and Xr — Vr are g au g e invariant, so it is 
physically meaningful to fix them to some value - this corresponds to locking the phase of 



the order parameter to those of the on-site s-wave density, i.e. the phase of the local Mott 
singlets. Similarly, by virtue of Eq. E3 the combination e % ^ Xr ~ XrtS} h r>r i is also gauge invariant. 
In the following, we shall set Xr + Xr' ~ Vr,r> — Xr — Vr = 0. 

The full set of variational parameters are now Xr, which characterize the transformation 
U , together with the parameters {«}, used to characterize U' = e zSu . In the subsequent 
mean field calculation, these will be n, A r y, A r and h r>r >, where A r is the amplitude of the 
local s-wave order parameter. 

It is now clear that in order to evaluate the expression in Eq. E3 we must have some 
information about the phases Xr- in the undoped Mott state, they only contribute to an 
overall phase, and are thus completely random. If this would be true also in the doped 
state, both terms in Eq. El would average to zero because of the phase factor e l ^ Xr ~ Xr '\ In 
the doped state, however, it is reasonable to assume that some short range correlation is 
generated among the phases Xr- This would correspond to having 

( e «"e- ix ' / )=/(r-f ? ). (19) 

Really, this should be shown by including the phases Xr in the variational calculation, but 
this is technically very hard to do, so we shall instead simply assume Eq. El to hold and set 
f(r — r') = 1 for nearest neighbors. Under this assumption, Eq. El simplifies to, 

( h rr ' ) = ~ h r y (1 — h) — y A r y n + ... — -^h r hh 2 , (20) 

where for future reference we have included the lowest higher order term which couples 
linearly to A r . We therefore find that the Hamiltonian expectation value to be minimized 
is given by the sum of Eq. [TT] and Eq. EDI subject to total particle number given by 
Eq. EJ This effective Hamiltonian looks very much like an ordinary BCS Hamiltonian, 
corresponding to Eq. ^ but with one dramatic difference, namely the presence of a term 
proportional to A rr > n, as well as a higher order term which couples linearly to the s-wave 
pairing operator. 

IV. MEAN FIELD ANALYSIS 

The mean field Hamiltonian can be analyzed by several equivalent methods. In the spirit 
of what was just developed, we could e.g. make a Bogoliubov-Valatin canonical transfor- 
mation and minimize the energy of the retransformed vacuum. This variational procedure 



would precisely correspond to the canonical transformation e lSu alluded to earlier. The 
method we actually use generalizes easily to finite temperatures and arbitrary large number 
of terms in the polynomial expansion of the mean field Hamiltonian. It uses that the density 
matrix p = e^^ T '^i~\ H ~^ N i) minimizes the free energy F = (H — pN) p — kT(S) p for all 
values of p (e _/3n is the partition function and S the entropy). We can then take p to be the 
exponential of an expression linear in n r , A rr /, A r and h r y , where the prefactors are varied 
to minimize F. This method yields the ordinary BCS theory when applied to a Hamiltonian 
of the form in Eq. ^ and gives a more complicated self-consistent calculation when more 
terms are kept. 

We have performed the mean field analysis numerically, both using the truncated expres- 
sions given explicitly above and the full mean field theory containing polynomials to seventh 
order. Since we construct an effective Hamiltonian, we define the hatted operators whose 
expectation values give the values AJ in Eq. O The corresponding operators A£, A* rr ,,h 
are therefore formally AJ = e~ tSu Aje* Gu etc. but in the calculation this involves simply 
reinterpreting the original operator in Eq. El in terms of quasiparticle fermion operators. 
For the density matrix p oc e~@ He ff we choose the Hamiltonian H = H e /f as 

Heff = J>' ^ + 2§ 7 Kr>) ~ 5>'^ + 2§ ^ A -) + ( CC ) ( 21 ) 

Our approximation is reasonable for small doping, and we confine the mean field analysis 
to this regime. We let e^ = (2 — cos k x — cos k y ) in two dimensions with a similar expression 
for D = 3. With p! = Dt' + 5 P and 7 2 = Dj + 5 7 we define e^ = t' e*. + 5^, d k = 7^ + 5 7 
and Ek = v/e| + d|. The energy gap A is the minimum in E^, and it is easy to verify 
that for small doping an excellent approximation is given by A = — p — = — . We define the 

momentum space sums as f[k) = ^2 k £ ■ Taking into account that there are four spin 
values we find the following expression for the doping S ~ n, where expectation values of 
operators are dropped when the context is clear, 

n = 4 x ^J (l - |J d 2 k = 2 -2e k /\t'\ (22) 

The expressions for if and 7 can be read off to lowest order from Eq. I2H 

t' = _i£(i_n) (23) 

5ht 

10 



7= ~m, (24) 



while the definition of the on-site s-wave order parameter can be read off from Eq. ^2 

L>7 - 5 7 = 2 v / 2f/A r . (25) 

In writing Eq. |^ we neglected the non-local repulsive interactions of type n r n r i that will 
certainly be present in any realistic model with screened Coulomb interaction. Assuming a 
nearest neighbor term, U\n r n r i and using the identity (n r n r >) = 1 + ^ A r , r 'A^ r/ — jh r yh r ' tr — 
2(n r + n r i) + |(nj! + n*,) + n r n r i, Eq. |2U would change to 

5ht 



7 



+ y/2U 1 K\ r , (26) 



4^2 

Below we argue that this would not qualitatively change our conclusions. 
As usual, self-consistency implies a gap equation which here reads, 

1 



-4 • (27) 



After expanding to lowest order in 5 and doing some some algebra this can be recast as 

J_ (A_ _ 5D6\ _ 55 _ _A~ 

which is a closed equation that can be used to find the physical gap A as a function of 8. 

Using the self-consistent equations, we can find expressions for the extended and onsite 
s-wave pairing amplitudes: 

' 5DS AA (29) 
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Thus, not surprisingly we find that the onsite s-wave component is reduced by a factor t/U 
relative to the extended component. We note that the particular combination corresponding 
to A rr > occurs in the left hand side of Eq. I2H1 which is consistent with the onsite s-wave 
component vanishing while the nearest neighbor extended pairing remains finite in the limit 
U — »• oo. We can now understand what would be the qualitative effect of adding extra 
repulsive interactions corresponding to the redefinition Eq. |2S1 of the variational parameter 
7. For large Lq, 7 will effectively be put to zero corresponding to A r>r r ~ ^n, rather than 
A r y ~ 5. We see that the scale of A r y changes but it is still non-zero for arbitrary small 
doping. 

11 



A. Asymptotic behavior of the gap in D = 2 

In two dimensions, the expression Eq. E21 can be approximated by 

n^^+y^). (30) 

The expression for 1 is logarithmically singular but can be approximated by 

T M zl h f -J±±VS±E\ = J_ ,„ CH^3!) , ( 3i) 

2vr I 32|f| J 2vr V A / 

where Eq. 0H1 was used. In spite of the logarithmic singularity, the self-consistent equation 
Eq. |28jcan be solved in closed form. Defining q = A/(t5) the (inverse) equation is 

87J- ir(4g-5 v / 2(l+2[f/t)) 

S&— e ^W* . (32) 

q 2 

By plotting the pairs (5, qS) according to the above formula as a function of q, we find the 

gap as a function of 5, shown for values of U — (oo, 10, 5, 1) in Fig. 1. We can see an almost 

linear behavior of the gap as a function of doping that is quite insensitive to the value of U. 

For all values of 5 and U the approximation A « t5 is a surprisingly good approximation. 

A comparision with the numerical solution of the self-consistent equations is shown in Fig. 

2. 

B. Asymptotic solution for D = 3 and small doping 

In the case D — 3, the vanishing density of states near k = makes the integral 1 
converge. In this case the self-consistent equation is Eq. EH] with D = 3, and constant 1. 
However, due to the vanishing density of states, even relatively small values of density lead 
to quite substantial values of jjl and 5 so the asymptotic value of this equation is far from 
being reached even for doping as low as .01. The corrections to 1 are rather slowly varying, 
so the linear dependence of the gap upon doping is obtained for the D = 3 case as well, as 
shown in Fig. 2. 

V. EFFECTIVE THEORY FOR SMALL 5 

After applying a transformation that rigorously preserves the full Fock space, we have 
obtained a non- Fermi liquid behavior for the doped spin 3/2 Mott insulator by using standard 

12 



mean field theory methods. Here we contrast this to an effective theory for the Hamiltonian 
in Eq. Q derived in the limit of small doping 5 <C 1 and small hopping t <C U, J, using 
the more conventional approach of projection on a low energy subspace. Ignoring hopping, 
t = 0, this low energy sector consists of states where each site is occupied by either a Mott 
singlet | $ s ), which we choose as the Fock vacuum, | $ s ) — > | ), or by a single charge 
| s ) = c\\ ) with spin 3/2 (s = ±1/2, ±3/2) - all other states are separated from these by 
a gap of order X ~ U, J. Restricting to this low energy sector and including the hopping 
in perturbation theory gives to order t 2 /X, the effective Hamiltonian 

H e ff = t 2 j C sr C s,r+8 + / J Jaa'(3[3'C ar C a i r Cp r+8 Cp^ r+ $ (33) 

subject to the constraint J2 s c tr c sr < 1 (Jaa'pp' are SU{2) scalars). This t — J-type model 
describes four species of fermions with nearest neighbor hopping (t ~ t), nearest neighbor 
exchange couplings (J ~ t 2 /X) and with the hard-core constraint that no two fermions 
occupy the same site. Second order perturbation theory guarantees a finite, albeit weak, 
attraction which opens the possibility of having a superconducting phase even for small 
doping. From this approach, however, we would expect such a phase to be destroyed by a 
nearest neighbor repulsion that is normally present in a realistic model. Thus, our previous 
mean field calculation is at odds with this approach. If the former turns out to be valid, 
it suggests that there are non-perturbative effects due to the hard-core constraints that are 
not easily accounted for in the conventional formulation. If on the other hand the hard-core 
constraints are not very important and the naive picture of four different species of weakly 
interacting fermions is essentially correct, it would suggest that our mean field treatment of 
the phase phase fluctuations does not capture the correct physics. 

VI. DISCUSSION AND SUMMARY 

A. The anomalous term A r yn r . 

We see from Eq. EI] that in order to have a nonzero doping 5 = 2— (n), we must have 
77- > 0, in fact n w 5 to lowest order in 5. Energetically we will also have h ^ for finite 
5. According to Eq. E3 the extended s-wave pairing field A r y cannot vanish and in fact 
will be proportional to doping. This in turn generates a (much smaller) on-site pairing 

13 



A,,, through the self-consistent equations. Note that this pairing field can never completely 
vanish because of the linear coupling to higher order terms. At finite temperature the 
mean field theory will presumably eventually break down via an xy-transition due to phase 
fluctuations that we have not taken into account. This has been discussed in a series of 



recent papers where the term "gossamer superconductivity" la, |l7 



Il8j ] has been used to 



describe a similar scenario. 

It is admittedly not easy to understand the physical origin of these new anomalous terms 
of the type A rr /n r . On a technical level, they are forced by the fermion statistics which 
constrains the form of the canonical transformation necessary to map a local Mott ground 
state to the vacuum. In our case, this transformation must be (a) nonlinear in fermion 
operators and (b) charge non- conserving. Property (a) yields an effective interaction from 
the hopping term near a charged ground state and property (b) makes this interaction non- 
gauge invariant. Property (a) is a necessary consequence of mapping a locally entangled state 
to the vacuum and property (b) is a conseqence of mapping a charged state to the vaccum 
which breaks gauge symmetry. Very general arguments relying on long range phase coherence 
and a finite range gap function then predict that the system should be a superconductor 
[11| . Our mean field calculation, which suggests a superconducting ground state, supports 
this picture, given our assumptions about phase coherence. 

We already pointed out the contradiction between our main result and what would be 
expected based on a conventional analysis of the type leading to Eq. ESI but also stressed 
the difficulties related to the hard core constraints inherent in the latter approach. Here we 
should note that more elaborate schemes for dealing with these non-holonomic hard core 
constraints face severe difficulties related to phase fluctuations. For example, in the spin 
1/2 Hubbard model at half-filling, one can turn the no double occupancy constraint into a 
holonomic gauge constraint by introducing spinons and holons. The resulting phase depends 
crucially on the fluctuations in the related gauge fields. By working in the full Hilbert 
space, we avoid these difficulties, but nevertheless our conclusions are still dependent on 
certain assumtions about phase coherence. Without a more sophisticated analysis of the 
phase fluctuations, we cannot rule out that these will be important and e.g. destroy the 
superconduting state at low doping. 
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B. Range of validity and applicability 

We now assume that our analysis is correct at low doping, and discuss its range of validity 
and applicability At sufficiently large value of doping, the theory will yield a free energy 
which is unfavorable compared to that of a doped n = 1 Mott insulator. The mean field 
picture suggests there will be a coexistence region where a slightly hole doped n = 2 Mott 
insulator will coexist with a hole doped n = 1 Mott insulator. The n = 1 Mott insulator 
will presumably have some sort of magnetic order at low temperature that breaks the large 
spin degeneracy of the uncorrelated odd filling Mott insulator. If a coexistence region really 
exists, or whether an intermediate phase which breaks translational invariance may exist, is 
beyond the scope of the present analysis. Our calculation thus makes assumptions about U 
which leave open the question if this behavior could really be seen in a physical system. On 
the one hand, U must be large enough (and J even larger) so that a Mott insulating state 
occurs at n = 2 and furthermore triply occupied sites are effectively absent. On the other 
hand, U must be small enough so that the correlated state will have lower energy than a 
mixed state with ann = l Mott insulator and ann = 2 Mott insulator. 

We should also ask how dependent our approximations are on our specific choice of 
parameters. In particular what would happen if we had taken J < in the spin 3/2 model? 
According to the analysis in Ref. [12j, the ground state at n = 2 is than 5*0(5) invariant 
and we could choose for instance c\ c\ |0) as the Mott vacuum. However, this as well as 
the other states with m ^ are not entangled, and the procedure discussed here results in 
a simple particle-hole transformation of pairs of fermion operators which does not generate 
any anomalous terms and the ordinary BCS-type analysis presented in these calculations 
should be valid. However, the state m = is entangled, and it can be checked by the 
methods of Ref. [15J that the canonical transformation that maps this state to the vacuum 
will be nonlinear in fermion operators as well as charge nonconserving causing an effective 
interaction and superconductivity to appear at finite doping through a similar mechanism 
discussed here for the Mott singlet. 

If we take the parameters as < J < 2U, which of course is more physically reasonable, 
the energetically most favorable state constructed with two quasicharge operators c\ r from 
the Mott singlet will no longer be a singlet with relative charge minus two, but rather one 
of the I = 2 multiplets with zero relative charge. The extended s-wave order parameter 
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will still appear since it is a direct consequence of mapping the Mott singlet to the vacuum, 
but there will be a spin order parameter introduced that takes the place of the local s-wave 
superconducting order parameter. The presence of an additional order parameter together 
with spin degeneracy makes this calculation more difficult and it awaits further analysis. 

Finally, it is relevant to ask whether the transformation used for the spin 3/2 case could 
be applied to the spin 1/2 systems. First consider the canonical transformation which maps 
between these the n = and n = 2 states. This is the ordinary particle-hole transformation 
which is not charge conserving. However, the doubly occupied singlet is created by c! r cj r |0) 
and hence is factorizable in fermion operators. The canonical transformation is therefore 
linear and no new interaction terms are introduced in the transformation. The physical 
properties of the system are symmetric under charge conjugation, which is sufficient for the 
particle-hole transformation not to generate any new behavior and the present analysis is 
uninteresting. In the case of the half-filled Hubbard model, the Mott ground state corre- 
sponds to one electron per site. This cannot be mapped to the vacuum through a canonical 
transformation without globally violating the Fermi anticommutation relations! 15]. 

C. Summary 

We have presented a new type of canonical transformation for the half-filled spin 3/2 
Hubbard model that maps the Mott insulator at half filling to the vacuum. This canonical 
transformation is straighforward to generalize to multi-band Hubbard models with a local 
spin singlet Mott insulating ground state. At finite doping, a self-consistent mean field 
theory for such a system results in a phase with long range phase coherence. An order 
parameter that is usually identified with extended s-wave superconductivity appears and is 
proportional to doping. The calculation appears to be in contradiction to other methods of 
attacking these kinds of problems, and we pointed out the difficulties with both approaches. 
Our calculations bear a striking resemblance to the "gossamer superconductor" scenario that 
has been recently introduced by Lauglin and coworkers. Although we have only explored a 
specific half-filled spin 3/2 Hubbard model, we believe that our method could be useful for 
variety of similar models with locally entangled Mott insulating ground states. 
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FIG. 1: Values of A/t as a function of doping 5 for U = (oo, 10, 5, 1) for 5 < 0.1, with the gap 
decreasing monotone with increasing U at given S. 
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FIG. 2: Values of A/t as a function of doping 5 for U/t = 10 for 5 < 0.05. Th upper curve is for 
D = 3 and the lower data points for D = 2. The solid line is the fit to the 2D asymptotic curve 
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